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Abstract 


This  study  focused  on  the  development  of  shear-lip  in  DT  specimens  under  quasi-static  loading 
rate.  The  objective  of  the  study  was  to  demonstrate  the  feasibility  of  a  numerical  simulation  of 
shear-lips  for  a  25mm  base  metal  DT  specimen  in  three-point  bend  test.  Three  real  stress-real 
strain  curves  representing  three  different  materials  with  varying  degrees  of  plastic  behaviour  were 
considered.  These  material  curves  are  designated  as  Curves  1,  2  and  3.  Curve  1  represents  the 
actual  real-stress-  real-strain  curve  of  the  350WT  steel  material.  On  the  other  hand,  Curve  2 
represents  the  same  curve,  but  with  the  fracture  strain  reduced  to  half  that  of  Curve  1 ;  while  Curve 
3  represents  a  similar  curve  as  Curve  1,  but  with  the  yield  stress  reduced  to  half  the  Curve  1  value. 
Crack  advance  and  element  or  node  release  was  based  on  two  possible  criteria,  namely  (a)  when 
the  average  strain  in  the  crack  tip  elements  that  occupy  a  physical  surface  of  1.25  mm  (0.050 
inches)  exceeds  the  strain  at  fracture  on  the  supplied  data;  and  (b)  when  the  strain  energy  density 
in  elements  occupying  a  volume  defined  by  1 .25  mm  (0.050  inches)  of  free  surface  at  the  crack  tip 
and  1.25  mm  (0.050  inches)  of  material  depth  exceeds  the  strain  energy  density  calculated  from 
the  supplied  real  stress  real  strain  curves. 

It  was  shown  that  the  plastic  strain  contour  patterns  obtained  from  the  finite  element  analyses 
could  be  used  to  estimate  the  shear-lip  sizes.  At  any  given  crack  size  (a/W  level),  the  shear-lip 
size  is  smallest  for  Curve  2  material  and  largest  for  Curve  3.  This  trend  is  reasonable  as  the  Curve 
2  material  is  the  least  plastic  (100%  strain  to  failure  as  opposed  to  200%  strain  for  Curves  1  and 
3);  and  the  Curve  3  material  is  the  most  plastic  (has  lower  yield  stress  of  197MPa  compared  to 
395MPa  for  Curves  1  and  2,  and  high  failure  strain  of  200%).  For  a  given  material,  the  shear-lip 
size  tends  to  reduce  with  increasing  crack  size.  For  instance,  for  the  Curve  1  material,  the  shear- 
lip  sizes  for  a/W  =  0.4, 0.45  and  0.5  are  1 0.95mm,  1 0.30mm  and  9.70mm,  respectively.  The  same 
trends  of  shear-lip  sizes  were  observed,  regardless  of  whether  the  plastic  strain  criterion  or  strain 
energy  density  criterion  was  used  for  crack  advance.  However,  for  any  given  material  and  crack 
size,  the  shear-lip  size  predicted  based  on  the  strain  energy  density  criterion  was  greater  than  that 
predicted  based  on  the  plastic  strain  criterion. 
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Resume 


Les  principaux  travaux  de  la  presente  etude  portaient  sur  la  croissance  du  debord  de  rupture  dans 
des  eprouvettes  de  torsion  ductile,  pour  une  vitesse  de  mise  en  charge  quasi  statique.  Ils  visaient  a 
demontrer  I’applicabilite  de  la  modelisation  numerique  des  debords  de  rupture  dans  le  cas  d’une 
eprouvette  de  torsion  ductile  en  metal,  de  25  mm,  soumise  a  un  essai  de  pliage  en  trois  points.  On 
a  examine  les  cas  des  courbes  effort  reel-deformation  reelle,  appelees  Courbes  1,  2  et  3, 
representant  trois  materiaux  caracterises  par  des  comportements  plastiques  distincts.  La  Courbe  1 
represente  la  courbe  effort  reel-deformation  reelle  experimentale  obtenue  pour  Lacier  350WT.  La 
Courbe  2,  quant  a  elle,  correspond  au  meme  type  de  courbe,  avec,  toutefois,  une  valeur  reduite  de 
deformation  a  la  rupture,  afm  d’obtenir  une  valeur  equivalente  a  celle  de  la  Courbe  1 .  Dans  le  cas 
de  la  Courbe  3,  qui  est  semblable  a  la  Courbe  1,  c’est  la  limite  apparente  d’elasticite  qui  a  etc 
reduite  de  moitie  par  rapport  a  celle  de  la  Courbe  1 .  La  propagation  de  la  fissure  et  la  relache  de  la 
contrainte  appliquee  aux  elements  ou  aux  noeuds  etaient  basees  sur  deux  criteres  possibles  :  a) 
lorsque  la  deformation  moyenne,  dans  les  elements  de  I’extremite  de  fissure  qui  occupent  une 
dimension  materielle  de  1 ,25  mm  (0,050  po),  est  superieure  a  la  deformation  a  la  rupture,  selon  les 
donnees  foumies,  etb)  lorsque  la  densite  d’energie  de  deformation,  dans  des  elements  occupant  un 
volume  delimite  par  une  dimension  libre  de  1 ,25  mm  (0,050  po)  a  Lextremite  de  fissure  et  par  une 
profondeur  de  1,25  mm  (0,050  po),  est  superieure  a  la  valeur  de  la  densite  d’energie  de 
deformation  calculee  a  partir  des  courbes  effort  reel-deformation  reelle  foumies. 

II  a  ete  demontre  que  les  courbes  modeles  de  deformation  plastique  obtenues  grace  a  Lanalyse  par 
elements  finis  peuvent  servir  a  estimer  la  taille  des  debords  de  mpture.  Pour  toute  taille  de  la 
fissure  (niveau  de  aAV),  la  taille  du  debord  de  rapture  est  minimale  dans  le  cas  du  materiau  de  la 
Courbe  2  et  maximale  dans  le  cas  de  celui  de  la  Courbe  3.  Cette  observation  est  acceptable,  car  le 
materiau  de  la  Courbe  2  est  le  moins  plastique  des  trois  (100  %  de  deformation  a  la  rupture,  par 
rapport  a  200  %  de  deformation  a  la  rupture,  dans  le  cas  des  Courbes  1  et  3),  et  le  materiau  de  la 
Courbe  3  est  le  plus  plastique  (il  est  caracterise  par  une  faible  limite  apparente  d’elasticite  de 
197  MPa,  comparativement  a  395  MPa,  pour  les  Courbes  1  et  2,  et  par  une  deformation  a  la 
rupture  elevee,  soit  200  %).  Pour  un  materiau  donne,  la  taille  du  debord  de  rupture  tend  a  diminuer 
en  function  de  I’accroissement  de  la  taille  de  la  fissure.  Ainsi,  dans  le  cas  du  materiau  de  la 
Courbe  1 ,  la  taille  des  debords  de  rupture  est  respectivement  de  1 0,95  mm,  1 0,30  mm  et  9,70  mm, 
pour  des  valeurs  de  aAV  de  0,4,  0,45  et  0,5.  Ces  tendances  relatives  a  la  variation  de  la  taille  des 
debords  de  rupture  ont  ete  observees  dans  le  cas  des  deux  criteres  susmentionnes  de  propagation 
de  la  fissure,  soit  celui  de  la  deformation  plastique  et  celui  de  la  densite  d’energie  de  deformation. 
Toutefois,  pour  tout  materiau  donne  et  toute  taille  donnee  de  la  fissure,  la  valeur  prevue  de  la  taille 
du  debord  de  fissure  basee  sur  le  critere  de  la  densite  d’energie  de  deformation  etait  toujours 
superieure  a  celle  prevue  pour  le  critere  de  la  deformation  plastique. 
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Executive  summary 


Introduction  or  Background 

Over  the  past  two  decades,  Defence  R&D  Canada  has  conducted  research  in  fracture 
mechanics  and  fracture  control  and  has  identified  shear  lip  as  the  most  economical  and 
practical  material  property  for  assessing  fracture  toughness  and  location  in  the  transition  curve 
for  steel  and  weldments. 

In  1991  DRDC  Atlantic  first  published  their  theory  that  the  plastic  zone  at  initiation  of  crack 
advance  and  the  developing  shear  lip  were  approximately  equal.  Over  the  past  ten  years 
DRDC  has  tried  to  verify  this  experimentally  and  in  the  past  two  years  has  supported  work  in 
several  locations  to  try  to  verify  this  numerically  with  finite  element  models. 

This  contract  was  Martec’s  first  attempt  to  determine  the  plastic  zone  shear  lip  relationship. 

Results 

Martec’s  model  was  not  fine  enough  to  provide  a  definitive  plastic  zone  size 
at  crack  initiation.  However,  comparison  with  their  work  on  shear  lip  and  stretch  zone  showed 
that  the  plastic  zone  size  obtained  for  the  25  mm  specimen  was  about  9  mm  at  DL/B=0.03; 
whereas,  from  the  present  study,  the  shear  lip  size  is  9.7-10.95  mm.  If  it  can  be  established  that 
crack  initiation  starts  around  DL.B=0.03,  then  their  analysis  supports  the  DRDC  hypothesis. 

Significance 

This  is  a  major  element  of  fracture  control  theory  and  its  verification  by  FE  greatly  simplifies 
the  approach  to  determining  risk  of  fracture  in  naval  vessels. 

Future  Plans 

To  help  this  approach  gain  wide  acceptance  DRDC  hopes  to  refine  the  approach  to  shear  lip 
determination  in  additional  configurations.  Accordingly  further  FE  model  work  is  planned. 
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Sommaire 


Introduction  on  historique 

Au  cours  des  deux  demieres  decennies,  R  &  D  pour  la  defense  Canada  (RDDC)  a  effectue  des 
travaux  de  recherehe  dans  les  domaines  de  la  mecanique  de  la  rupture  et  sa  prevention.  Les 
resultats  ont  permis  d’etablir  que  pour  evaluer  la  resistance  a  la  rupture  (ou  tenacite)  et  la 
position  de  la  fracture  dans  la  courbe  de  transition  de  I’acier  et  des  assemblages  soudes,  le 
debord  de  rupture  constitue  la  caracteristique  materielle  dont  1’ analyse  est  la  plus  economique 
et  pratique. 

En  1991,  RDDC  Atlantique  publiait  sa  theorie  selon  laquelle  la  taille  de  la  zone  plastique  a 
I’amorqage  de  la  propagation  de  la  fissure  et  celle  du  debord  de  rupture  en  developpement 
etaient  approximativement  egales.  Depuis  dix  ans,  RDDC  tente  d’apporter  une  confirmation 
experimentale  a  la  theorie  et,  au  cours  des  deux  demieres  annees,  elle  a  soutenu  les  travaux  de 
differents  centres  de  recherche  visant  la  validation  numerique  de  la  theorie  a  partir  de  modeles 
d’elements  finis. 

Le  present  contrat  constituait  les  premiers  travaux  de  la  societe  Martec  visant  a  determiner  la 
relation  entre  la  zone  plastique  et  le  debord  de  rapture. 

Resultats 

Le  modele  utilise  par  la  societe  Martec  n’etait  pas  assez  raffine  pour  permettre  de  bien 
determiner  la  taille  de  la  zone  plastique  a  I’amor^age  de  la  fissure.  Toutefois,  en  comparant  les 
resultats  avec  ceux  obtenus  lors  de  leurs  travaux  portant  sur  le  debord  de  rupture  et  la  zone 
d’etirement,  on  a  etabli  que  la  taille  de  la  zone  plastique  etait  d’ environ  9  mm,  dans  le  cas  de 
I’echantillon  de  25  mm,  pour  un  rapport  DL/B  =  0,03.  Cependant,  dans  la  presente  etude,  la 
taille  du  debord  de  rupture  se  situe  entre  9,7  et  10,95  mm.  Si  les  resultats  demontrent  que 
I’amorqage  de  la  fissure  correspond  a  une  valeur  approximative  de  DL/B  =  0,03, 1’analyse 
soutient  I’hypothese  de  la  RDDC. 

Importance  des  resultats 

II  s’agit  d’un  element  majeur  de  la  theorie  de  la  prevention  des  ruptures.  Sa  verification  par 
1’analyse  par  elements  finis  simplifiera  grandement  les  methodes  de  determination  des  risques 
de  rupture  dans  les  navires  de  guerre. 

Futurs  travaux 

Pour  que  cette  methode  soit  largement  acceptee,  RDDC  souhaite  la  raffmer  et  I’utiliser  pour 
determiner  les  debords  de  rupture  dans  des  configurations  supplementaires.  Par  consequent,  on 
prevoit  executer  d’autres  travaux  de  modelisation  par  elements  finis. 
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1. 


INTRODUCTION 


1. 1  Background 

A  major  component  of  the  DRDC  Atlantic  Dockyard  Laboratory  (DRDC/ Atlantic  DL) 
fracture  control  program  is  linking  large  specimen  size  dynamic  tear  (DT)  transition  curves  to 
structural  transition  curves  and  quantitative  fracture  parameters.  To  do  this  it  is  necessary  to 
determine  relationships  between  stretch  zone  and  plastic  zone  size  before  crack  advance  and 
plastic  zone  size  at  and  just  after  crack  advance  and  the  size  of  the  developing  shear-lip.  The 
former  relationship,  that  between  stretch  zone  and  plastic  zone,  has  been  investigated  in  previous 
studies  resulting  in  a  relationship  of  plastic  zone  to  stretch  zone  of  80  to  100.  In  several  papers 
going  back  to  1991,  DRDC  Atlantic  DL  has  estimated  that  the  shear-lip  (actual  shear- lip  width 
plus  specimen  contraction)  that  develops  is  approximately  equal  to  the  mid  specimen  plastic  zone 
at  crack  advance. 

1.2  Objectives  and  Scope 

The  objective  of  this  study  is  to  demonstrate  the  feasibility  of  a  numerical  simulation  of 
shear-lips  for  a  base  metal  quasi-static  loading  rate  dynamic  tear  (DT)  specimen  three-point  bend 
test.  Three  real  stress-real  strain  curves  acquired  for  gauge  lengths  of  1.25  mm  (0.050  inches) 
supplied  by  DRDC/DL  will  be  used  for  the  model.  Crack  advance  and  element  or  node  release 
will  be  based  on  two  possible  criteria.  The  two  criteria  are  (a)  when  the  average  strain  in  the  crack 
tip  elements  that  occupy  a  physical  surface  of  1.25  mm  (0.050  inches)  exceeds  the  strain  at 
fracture  on  the  supplied  data;  and  (b)  when  the  strain  energy  density  in  elements  occupying  a 
volume  defined  by  1.25  mm  (0.050  inches)  of  free  surface  at  the  crack  tip  and  1.25  mm  (0.050 
inches)  of  material  depth  exceeds  the  strain  energy  density  calculated  from  the  supplied  real  stress 
real  strain  curves. 
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2. 


PROBLEM  DESCRIPTION 


2.1  Specimen  Configurations 

A  schematic  representation  of  the  test  specimen  is  presented  in  Figure  2.1.  The 
specimen  dimensions  are  presented  in  Table  2.1.  The  specimen  is  made  of  350WT  steel  material. 

2.2  Description  of  Shear-Lip  and  Plastic  Zone  in  DT  Specimens 

Figure  2.3  shows  schematically  the  plastic  zone  associated  with  the  DT  specimen.  In  this 
figure  Szw  is  the  stretch  zone  width  and  Vy  is  the  plastic  radius.  Figure  2.4  is  a  schematic  of  the 
crack  front  of  a  standard  DT  specimen  showing  the  shape  of  the  shear-lip,  S.  The  shear  index  can 
be  calculated  using  Equation  2.1,  where  SI  is  the  shear  index,  S  is  the  shear-lip  and  B  is  the 
specimen  width  (25mm  for  this  study).  [4] 

2S 

SI  =  —  Equation  2. 1 

2.3  Requirements 

It  is  required  to  model  the  shear-lip  formation  and  plastic  zone  development  by  finite 
element  technology  for  25  mm  standard  DT  specimens  with  cracks  extended  to  14  mm  by  fatigue 
as  the  starting  point.  Three  different  1 .25mm  gage  length  real  stress  real  strain  curves  will  be  used 
to  model  the  material  behaviour  and  define  element  or  node  release.  Formation  of  the  crack  and 
subsequent  shear-lips  at  a/W  of  0.4,  0.5,  and  0.6,  will  be  compared  to  the  plastic  zone  at  crack 
initiation  and  at  .050  inch  crack  extension.  The  specific  goals  of  the  study  include: 

1 .  Analysis  of  the  three  real  stress-real  strain  curves  for  the  data  supplied  to  determine  the 
strain  energy  density  for  each. 

2.  Review  of  three  available  approaches  for  modeling  the  crack  advance  -  element  removal, 
remeshing  and  node  release,  to  determine  which  one  is  most  appropriate  to  use  in 
conjunction  with  the  FE  program  chosen. 

3.  Development  of  FE  models  for  the  25  mm  DT  specimen  configuration,  model  element 
release  based  on  strain  to  fracture  and  strain  energy  density  and  determine  the  shear-lip 
formation  for  each  data  set  and  each  failure  criteria.  These  results  will  be  used  to 
demonstrate  the  sensitivity  of  the  simulated  shear-lip  size  to  the  modeling  parameters  and 
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failure  criteria. 

4.  Provide  results  in  tabular  and/or  graphical  form  to  display  the  essential  relationship  of 
plastic  zone  at  initiation  of  crack  advance  and  plastic  zone  at  1 .25mm  crack  advance  and 
the  developing  shear-lip  at  a/W  of  0.4  0.5  and  0.6. 

5.  Provide  motion  clips  for  viewing  dynamically  these  results. 

6.  Documentation  of  results  in  a  contractor  report. 

The  material  test  data  to  be  used  (1.25mm  gage  length  real  stress-real  strain  data)  is 
presented  in  Table  2.2  and  Figure  2.2.  Curve  1  represents  the  complete  real  stress-real  strain 
behaviour,  with  the  failure  strain  at  206%  and  yield  stress  at  394MPa.  Curve  2  is  similar  to  Curve 
1  except  that  the  failure  strain  is  103%,  which  is  half  of  the  Curve  1  value.  Similarly,  Curve  3  is 
the  case  with  the  failure  strain  the  same  as  Curve  1  but  with  the  yield  stress  equal  to  197MPa, 
which  is  half  of  the  Curve  1  value.  As  stated  before,  these  will  be  used  to  study  the  sensitivity  of 
the  shear-lip  size  to  the  modelling  parameters  and  failure  criteria. 
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Table  2.1:  Standard  DT  Specimen  Configuration  Dimensions 


Dimensions 

Specimen  (DT-350WT-25) 

Beam  Length,  L  (mm) 

181 

Beam  Span,  S  (mm) 

164 

Crack  Length,  a  (mm) 

14* 

Beam  Width,  W  (mm) 

41 

Beam  Thickness,  B  (mm) 

25 

*a  =  12  mm  plus  0.254  mm  pressing  followed  by  fatigue  extension  to  14  mm 


Table  2.2:  Material  Test  Data  to  be  Used  in  Analyses 


Strain 

(mm/mm) 

Stress  (MPa) 

Curve  1 

Curve  2 

Curve  3 

0 

395 

395 

197 

0.0133 

509 

509 

255 

0.3330 

650 

650 

325 

1.0000 

820 

820 

410 

1.0300 

820 

374 

410 

1.0310 

820 

2.4 

410 

1.6600 

949 

475 

2.0000 

949 

475 

2.0600 

374 

187 

2.0610 

2.4 

1.2 

Figure  2.1:  Configuration  of  DT  Specimen 
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Real  Stress  Versus  Real  Strain 


Real  Strain  (mm/mm) 

— ♦ —  Curve  1  — ■ —  Curve  2  a  Curve  3 


Figure  2.2:  Graphical  Representation  of  Material  Test  Data  to  be  Used  in  Analyses 


Figure  2.3:  Schematic  of  DT  Specimen  Showing  Plastic  Zone 
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Figure  2.4:  Schematic  of  Cross-Section  of  DT  Specimens  Showing  Shear-Lip  Formation 
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3. 


FINITE  ELEMENT  MODEL 


3.1  Finite  Element  Approach 

Several  different  methods  of  simulating  eraek  advanee  in  a  finite  element  model  have 
been  investigated,  including  remeshing,  node  release  and  element  removal  [5],  Consideration  was 
given  to  both  ABAQUS  [6]  and  ANSYS  [7]  and  finite  element  packages  to  perform  the  analysis 
needed  for  this  study.  Due  to  the  complexity  and  the  computation  effort  required,  it  was 
concluded  that  the  method  of  remeshing  would  not  be  considered  in  this  study.  Initially,  both 
element  removal  and  node  release  were  considered  in  the  ANSYS  environment.  Focus  was  placed 
on  the  method  of  node  release.  This  process  required  a  fairly  high  degree  of  computational  effort 
due  to  the  requirement  of  fairly  small  mesh  sizes.  If  the  mesh  were  too  coarse,  the  solution  would 
begin  to  diverge. 

Due  to  difficulties  experienced  using  ANSYS,  focus  was  shifted  to  the  ABAQUS.  The 
method  of  element  removal  using  the  MODEL  CHANGE  option  in  ABAQUS  was  considered. 
The  software  also  has  a  wide  range  of  non-linear  materials  models  and  can  account  for  large 
strains  and  displacements.  Initial  testing  of  this  software  for  this  procedure  was  very  promising, 
therefore,  the  ABAQUS  finite  element  software  was  used  to  model  the  elastic-plastic  behaviour  of 
the  Standard  DT  specimen.  However,  for  convenience,  the  HyperMesh  [4]  general-purpose  pre- 
and  post-processing  program  was  used  for  model  generation  and  results  processing. 

3.1.1  Element  Removal  Process  in  ABAQUS 

Figure  3.1  shows  a  flow  chart  outlining  the  procedure  adopted  to  perform  the  crack 
propagation  analyses.  An  external  program,  written  in  the  Fortran  language,  is  used  to  automate 
iterative  removal  process.  This  program  starts  the  ABAQUS  runs  and  then  reads  the  ABAQUS 
ASCII  result  files  once  an  analysis  step  is  complete  and  determines  if  any  nodes  or  elements  on 
the  crack  front  have  exceeded  the  specified  criteria  for  crack  growth  (plastic  strain  or  strain  energy 
density).  If  any  nodes  have  failed  the  criteria,  the  program  generates  and  runs  a  removal  step 
containing  any  elements  in  the  crack  front  connected  to  those  nodes.  The  ASCII  results  file  is 
read  after  this  removal  step  to  determine  if  any  nodes  or  elements  still  exceeded  the  crack  growth 
criteria.  If  there  are  additional  elements  to  be  removed  another  removal  step  is  performed. 
Otherwise  the  load  is  incremented  and  the  removal  process  is  repeated.  This  process  is  continued 
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until  all  of  the  elements  in  the  crack  front  have  been  removed,  until  a  specified  maximum  load  has 
been  attained  or  a  specified  maximum  number  of  steps  has  been  performed. 

3.2  Finite  Element  Meshes 

The  HyperMesh  software  was  used  to  generate  the  finite  element  models  of  the  DT 
standard  specimens,  which  were  then  translated  to  ABAQUS  input  files.  Several  finite  element 
models  were  generated  to  determine  a  mesh  refinement,  which  minimized  computational  effort  but 
did  not  diminish  quality  of  the  results.  The  model  selected  is  shown  in  Figure  3 .2(a)  and  (b).  This 
is  a  doubly  symmetric  model  (about  longitudinal  and  transverse  axis)  and  contains  4450  8-node 
brick  elements  and  5446  nodes.  The  elements  to  be  removed  in  the  crack  front  are 
0.683x0. 625x0. 25mm,  which  will  allow  the  crack  to  propagate  0.683mm  (0.027in)  at  a  time. 

It  should  be  noted  that  the  problem  being  solved  is  a  highly  non-linear  one  that  required 
enormous  computational  time  and  data  storage  requirements.  For  these  reasons,  the  mesh  size  and 
configuration  had  to  be  limited  to  the  one  shown  in  Figure  3.2.  With  this  model,  a  typical  run 
requires  2-3  days  and  1 .0-2.0  GB  of  data  storage.  To  further  reduce  the  storage  requirements,  only 
stress/ strain  results  on  only  the  elements  in  the  crack  face. 

3.3  Material  Models 

An  incremental  rate  independent  plasticity  theory  available  in  the  ABAQUS  finite 
element  program  [3]  was  used  for  the  material  constitutive  model.  This  standard  model  for 
plasticity  is  summarized  in  the  Phase  I  report  [1].  As  previously  mentioned,  three  different 
representations  of  stress-strain  behaviour  of  350  WT  steel  material  are  considered.  These  curves 
were  presented  in  Table  2.2  and  Figure  2.2. 
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3.4  Boundary  Conditions  and  Loading 

The  following  boundary  conditions  were  applied: 

At  the  support:  y  =  0 

Along  the  longitudinal  centre  line:  z  =  0 

Along  the  transverse  centre  line:  x  =  0 ; 

where,  x,y  and  z  are  the  displacements  components  in  the  longitudinal,  transverse  and  out  of  plane 
directions,  respectively.  The  boundary  conditions  are  displayed  in  Figure  3.3. 

The  load  is  applied  as  a  pressure  load  distributed  across  a  8mm  section  at  the  top  of  the 
coupon  (darkest  elements  in  Figure  3.2)  as  shown  in  Figure  3.3. 
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Figure  3.1:  Flow  Chart  Showing  the  Procedure  Used  for  Element  Removal  in  ABAQUS 

Environment 
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(a) 


Crack  Front 
(a  =  14mm) 


(b) 


Figure  3.2:  Finite  Element  Mesh  of  Standard  DT  Speeimen 
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Longitudinal 

Symmetry 


Applied 
Pressure  Load 

Transverse 

Symmetry 


Simple  Support 


Figure  3.3:  Finite  Element  Model  of  Standard  DT  Speeimen  Showing  Applied  Loading  and 

Boundary  Conditions 
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4. 


RESULTS  AND  DISCUSSIONS 


4.1  Strain  Energy  Densities  for  the  Three  Curves 

The  three  real  stress  real  strain  curves  provided  were  analysed  and  the  corresponding 
strain  energy  densities  were  determined  as  follows  (by  calculating  the  area  under  the  stress  vs 
strain  curves): 

Curve  1  -  1625  N  mm/mm^ 

Curve  2  -  700  N  mm/mm^ 

Curve  3  -  813  N  mm/mm^ 

These  results  are  used  as  the  criteria  value  for  element  removal  using  strain  energy  density. 

4.2  Shear-Lip  and  Plastic  Zone  Development  with  Maximum  Strain  Criteria 

During  the  analysis  process  a  computer  disk  space  problem  was  encountered.  When 
storing  results  for  all  of  the  elements  in  the  model  for  each  analysis  step,  the  result  file  sizes  were 
quickly  exceeding  1 .0GB.  It  was  determined  that  there  was  no  significant  advantage  to  retaining 
results  for  any  elements  outside  of  those  in  the  plane  of  the  crack  front.  Figure  4.1  shows  the 
plastic  strain  contours  at  a  displacement  level  of  Dl  =  6.29mm.  As  can  be  seen  in  this  figure,  the 
plastic  zone  resulting  from  the  crack  front  can  no  longer  be  distinguished.  This  plastic  zone 
merges  with  the  yielded  region  from  the  load  after  a  displacement  of  approximately  1mm.  For  this 
reason,  results  are  only  printed  for  elements  in  the  plane  of  the  crack  front  (darker  elements  in 
Figure  4.2). 

Figure  4.3  shows  the  process  used  to  determine  the  shear-lip  size,  S  from  the  finite  element 
analysis.  As  illustrated  in  the  figure,  the  shear-lip  occurs  as  a  result  of  the  so-called  tunnelling 
effect  where  the  crack  grows  faster  in  the  centre  of  the  specimen  as  a  result  of  higher  stress 
triaxialilty  in  the  centre  than  at  the  edges.  Thus  the  material  tends  to  break  at  the  centre  while  the 
edges  try  to  hold  on.  The  shear-lip  size,  S  is  given  by  the  largest  distance  from  the  edge  of  the 
specimen  to  the  edge  of  the  fractured  surface  (“tunnel”)  in  the  centre  of  the  specimen.  To  obtain 
the  shear-lip  size  from  the  finite  element  analysis,  plastic  strain  contours  of  the  crack  face  are 
plotted  for  each  case.  The  fractured  surface  corresponds  to  the  region  with  very  high  plastic 
strains  shown  as  the  red  contours  in  the  figures  that  follow.  The  shear-lip  is  measured  from  the 
specimen  edge  to  the  point  where  the  high  strain  contours  plateau  out  near  the  specimen  centre,  as 
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shown.  This  process  was  shown  to  provide  the  most  consistent  and  reasonable  set  of  results. 
However,  in  some  cases,  it  was  difficult  to  measure  this  distance  due  to  big  jumps  in  contours 
obtained,  and  judgement  had  to  be  employed  in  determining  the  shear-lip  size.  These  cases  are 
highlighted  in  the  discussions  below. 

Figure  4.4  shows  the  plastic  strain  contours  and  the  corresponding  location  of  the  crack 
front  (degree  of  element  removal)  for  the  case  using  Curve  1  material  definition  and  a  plastic 
strain  of  200%  as  the  removal  criterion.  Figure  4.4(a)  shows  the  plastic  strain  contours  on  the 
crack  face  and  Figure  4.4(b)  shows  the  elements  removed  when  the  crack  front  has  reached  an 
a/W  value  of  0.4.  Figure  4.4(c)  shows  a  close-up  view  of  the  plastic  strain  contours  that  are  used 
to  estimate  the  shear-lip  size.  Figures  4.5  and  4.6  show  the  corresponding  plots  for  the  cases  when 
the  crack  front  level  reached  a/W  values  of  0.45  and  0.5,  respectively.  Unfortunately,  due  to 
convergence  problems,  the  analysis  could  not  be  continued  to  a/W  =  0.6. 

Figure  4.7  to  4.9  show  the  results  for  the  case  using  Curve  2  material  definition  and  a 
plastic  strain  of  100%  as  the  removal  criterion.  Figure  4.10  to  4.12  show  the  results  for  the  case 
using  Curve  3  material  definition  and  a  plastic  strain  of  200%  as  the  removal  criterion. 

The  shear-lip  size(s)  obtained  for  the  three  material  curves  are  summarized  in  Table  4.1. 
The  corresponding  shear  indices  (SI  =  2S/B)  are  also  presented  in  the  table,  where  B  is  the 
specimen  thickness.  At  any  given  crack  size  (a/W  level),  the  shear-lip  size  is  smallest  for  Curve  2 
material  and  largest  for  Curve  3.  This  trend  seems  reasonable  as  the  Curve  2  material  is  the  most 
brittle  (100%  strain  to  failure  as  opposed  to  200%  strain  for  Curves  1  and  3);  and  the  Curve  3 
material  is  the  most  plastic  (has  lower  yield  stress  of  197MPa  compared  to  395MPa  for  Curves  1 
and  2,  and  high  failure  strain  of 200%).  For  a  given  material  (Curve  1  for  example),  the  shear-lip 
size  tends  to  reduce  with  increasing  crack  size.  For  instance,  for  the  Curve  1  material,  the  shear- 
lip  sizes  for  a/W  =  0.4,  0.45  and  0.5  are  10.95mm,  10.30mm  and  9.70mm,  respectively.  This 
observation  is  consistent  with  earlier  DRDC  -  Atlantic  studies  [4]  in  which  it  was  shown  that  the 
shear-lip  size  and  hence  shear  index  (SI)  was  related  to  the  fracture  energy.  The  specimens  with 
larger  size  cracks  required  smaller  fracture  energy  and  hence  smaller  shear-lip  sizes.  A  similar 
trend  is  generally  observed  for  the  Curves  2  and  3  materials,  although  some  discrepancies 
occurred  at  a/W  =  0.4  for  Curve  2  and  at  a/W  =  0.45  for  Curve  3.  These  are  some  of  the  cases  for 
which  it  was  difficult  to  measure  the  shear-lip  sizes  precisely,  as  discussed  above. 
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4.3  Shear-Lip  Development  with  Strain  Energy  Density  Criteria 

The  use  of  the  strain  energy  density  criterion  presented  some  difficulties  with  regards  to 
convergence  of  the  solution.  This  was  probably  because  of  the  fact  that  the  fracture  criterion  was 
based  on  the  uniaxial,  rather  than  the  triaxial  strain  energy  density;  and  also  due  to  the  rather  large 
element  sizes  (about  0. 683x0. 625x0.25mm)  in  the  fracture  face.  Fracture  usually  occurred  at 
higher  strain  levels  than  the  fracture  strains,  hence  the  stress-strain  curve  had  to  be  extended 
somewhat  beyond  the  fracture  strain,  to  achieve  convergence.  Due  to  these  difficulties  as  well  as 
time  and  budget  constraints  it  was  not  possible  to  complete  all  of  the  cases  planned.  Only  analysis 
for  Curves  1  and  2  were  performed. 

Figures  4.13  to  4.14  show  the  plastic  strain  contours,  the  corresponding  location  of  the 
crack  front  and  close-up  views  of  the  strain  contours  for  the  case  using  Curve  1  material 
definition,  using  a  strain  energy  density  of  1625  N  mm/mm^  as  the  removal  criterion,  when  the 
crack  front  has  reached  a/W  or  0.4  and  0.45,  respectively.  Figures  4.15  to  4.17  show  the 
corresponding  cases  for  the  Curve  2  material,  using  a  strain  energy  density  of  700  N  mm/mm^  as 
the  removal  criterion,  when  the  crack  front  has  reached  a/W  of  0.4,  0.45  and  0.5,  respectively. 

The  shear-lip  sizes  obtained  for  these  cases  are  also  summarized  in  Table  4. 1 ,  where  the 
corresponding  shear  indices  are  also  presented.  Again,  the  shear-lip  sizes  are  smaller  for  the  more 
brittle  material  (Curve  2)  than  the  more  plastic  material  (Curve  1).  For  instance,  at  a/W  =  0.4,  the 
shear-lip  sizes  for  the  Curve  2  and  Curve  1  materials  are  1 1.90mm  and  8.10mm,  respectively. 
Furthermore,  for  a  given  material  the  shear-lip  size  decreases  with  increase  in  crack  size. 
Consider,  for  instance,  the  case  of  the  Curve  1  material.  The  shear-lip  size  varied  from  1 1 .90mm 
at  a/W  =  0.4  to  1 1 .20mm  at  a/W  =  0.45.  However,  this  trend  is  not  very  precise  with  the  Curve  2 
material. 

From  Table  4.1,  it  is  seen  that  for  any  given  material  and  crack  size,  the  shear-lip  size 
predicted  based  on  the  strain  energy  density  criterion  is  greater  than  the  predicted  based  on  the 
plastic  strain  criterion.  As  discussed  above,  this  is  probably  due  to  the  strain-energy  density  being 
obtained  from  uniaxial  testing  rather  than  a  triaxial  test,  such  that  fracture  occurred  at  higher  strain 
levels  than  the  fracture  strain.  A  more  refined  approach  will  be  to  use  strain  energy  density 
obtained  from  material  tested  in  a  triaxial  state  of  stress  that  is  closer  to  the  state  of  stress  in  the 
DT  specimen. 
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Table  4.1:  Summary  of  FE  Results  -  Shear-Lip  Sizes  for  25  mm  Thick  Specimen 


a/w 

Property 

Curve  1 

Curve  2 

Curve  3 

Strain 

Criterion 

Strain 

Energy 

Criterion 

Strain 

Criterion 

Strain 

Energy 

Criterion 

Strain 

Criterion 

Strain 

Energy 

Criterion 

0.40 

s  (mm) 

10.95 

11.90 

7.50 

8.10 

11.55 

- 

SI 

0.876 

0.952 

0.600 

0.648 

0.924 

- 

0.45 

s  (mm) 

10.30 

11.20 

8.70 

8.60 

10.00 

- 

SI 

0.824 

0.896 

0.696 

0.688 

0.800 

- 

0.50 

s  (mm) 

9.70 

- 

6.85 

8.60 

10.9 

- 

SI 

0.776 

- 

0.548 

0.688 

0.872 

- 

slepi  incIO,  l=G.09e-001 
Plastic  Strain  Equivalent 


>  l  .OOe-03 

<  1.00e-03 

<  9.23e-04 

<  8.46e-04 

<  7.G9e-04 

<  6.92e-04 

<  6.15e-04 

<  5.38e-04 

<  4.62e-04 

<  3.85e-04 

<  3.08e-04 

<  2.31e-04 

<  1.54e-04 

<  7.B9e-05 

<  O.OOe  +  00 


Max  =  7.30e-01 
Min  =  O.OOe  +  00 


Figure  4.1:  Plastic  Strain  Contours  at  Dl  =  6.29mm 
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Figure  4.2:  Finite  Element  Mesh  Showing  Elements  For  Which  Results  are  Stored 


Figure  4.3:  Schematic  of  Crack  Face  Showing  Method  of  Determining  Shear-Lip  Size 
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step30  inci,  l=2.13e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1 .50e  +  00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-81 

<  1.15e-B1 

<  8.88e  +  80 


Max  =  1.8Be  +  B0 
Min  =  8.80e  +  08 


A 


(a) 


(b) 


step30  inci,  1=2.1 3e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1.04e+00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-B1 

<  5.77e-81 

<  4.626-81 

<  3.46e-B1 

<  2.316-81 

<  1.156-81 

<  0.086  +  00 


Max  =  1 .806  +  00 
Min  =  0.006  +  00 


A 


<1- 

z 


(c) 


Figure  4.4:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  1  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.4) 
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step42  inc1,  t=2.94e-i-001 
Plastic  Strain  Equivalent 


>  1.50e  +  00 

<  1.50e  +  00 

<  1.38e  +  00 

<  1.27e  +  00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-01 

<  1.15e-01 

<  O.OOe  +  00 


Max  =  1.85e  +  BB 
Min  =  O.OOe  +  00 


A 


-<1- 

z 


(a) 


(b) 


step42  inci,  t=2.94e+001 
Plastic  Strain  Equivalent 


S>  1.50e+00 
<  1.50e+00 
<  1 .38e  +  00 
<  1.27e+00 
<  1.15e  +  00 
<  1 .04e  +  00 
<  9.23e-01 
<  8.08e-01 
<  6.92e-01 
<  5.77e-01 
<  4.62e-01 
<  3.46e-01 
<  2.31e-01 
<  1.15e-01 
<  O.OOe  +  00 

Max  =  1.85e  +  00 
Min  =  O.OOe  +  00 


A 


<3- 

z 


(C) 

Figure  4.5:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  1  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.45) 


DRDC  Atlantic  CR  2002-178 


19 


step58  inci,  l=3.95e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1  O^e+OO 

<  9.23e-01 

<  8.08e-01 

<  6.92e-B1 

<  5.77e-B1 

<  ^.62e-B1 

<  3.-16e-B1 

<  2.31e-B1 

<  1.15e-B1 

<  8.88e  +  88 


Max  =  1.89e  +  BB 
Min  =  8.88e  +  88 


A 


z 


(a) 


(b) 


step58  inci,  l=3.95e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1 .50e  +  00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.466-01 

<  2.31  e-01 

<  1  15e-01 

<  O.OOe+00 


Max  =  1.89e  +  00 
Min  =  O.OOe  +  00 


A 


<1- 

z 


(C) 

Figure  4.6:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  1  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.5) 
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stepIO  inc4,  l=7.65e+000 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  B.OOe-OI 

<  4.62e-01 

<  4.23e-01 

<  3.85e-01 

<  3.46e-01 

<  3.08e-01 

<  2.70e-01 

<  2.31 e-01 

<  1.93e-01 

<  1.54e-01 

<  1.16e-01 

<  7.75e-02 

<  3.91  e-02 

<  6.39e-04 


Max  =  1.59e  +  00 
Min  =  6.39e-64 


A 


O- 

z 


(a) 


(b) 


stepIO  inc4,  t=7.65e+000 
Plastic  Strain  Equivalent  (N) 


S>  5.00e-01 
<  5.00e-01 
<  4.62e-01 
<  4.23e-01 
<  3.85e-01 
<  3.46e-01 
<  3.08e-01 
<  2.7Be-01 
<  2.31e-01 
<  1.93e-01 
<  1.54e-01 

_  <1.16e-01 

_  <  7.75e-02 

_  <3.91e-02 

<  6.39e-04 

Max  =  1.59e  +  0B 
Min  =  6.39e-04 


A 


<]- 

z 


(C) 

Figure  4.7:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.4) 
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step14  inc4,  l=1.17e+001 
Plastic  Strain  Equivalent  (N) 


B>  5.00e-01 
<  5.00e-01 
<  4.62e-01 
<  4.24e-01 
<  3.86e-01 


<  3.10e-01 

<  2.72e-01 

<  2.34e-01 

<  1.96e-01 

<  1.58e-01 

<  1.20e-01 

<  8.18e-02 

<  4.38e-02 

<  5.77e-03 

Max  =  1.B3e  +  0B 
Min  =  5.77e-B3 


A 


<3- 

z 


(a) 


(b) 


step14  inc4,  t=1.17e+001 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  5.00e-01 

<  4.62e-01 

<  4.24e-01 

<  3.86e-01 

<  3.48e-01 

<  3.10e-01 

<  2.72e-01 

<  2.34e-01 

<  1.96e-01 

<  1.58e-01 

<  1.20e-01 

<  8.18e-02 

<  4.38e-02 

<  5.77e-03 

Max  =  1.G3e  +  00 
Min  =  5.77e-03 


A 


<]- 

z 


(c) 


Figure  4.8:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.45) 
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step18  inci,  l=1.57e+001 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  5.00e-01 

<  4.62e-01 

<  4.25e-01 

<  3.87e-01 

<  3.50e-01 

<  3.12e-01 

<  2.74e-01 

<  2.37e-01 

<  1.99e-01 

<  1.62e-01 

<  1.24e-01 

<  8.63e-02 

<  4.87e-82 

<  1.11e-82 


Max  =  1.35e  +  00 
Min  =  1.11e-02 


A 


O- 

z 


(a) 


(b) 


stepIS  inci,  t=1.57e+001 
Plastic  Strain  Equivalent  (N) 


>  5.00e-01 

<  B.OOe-OI 

<  4.62e-01 

<  4.25e-01 

<  3.87e-01 

<  3.50e-01 

<  3.12e-01 

<  2.74e-01 

<  2.37e-01 

<  1.99e-01 

<  1.62e-01 

<  1.24e-01 

<  8.63e-02 

<  4.87e-02 

<  1.11e-02 


Max  =  1.35e  +  00 
Min  =  1.11e-02 


A 


<3- 

z 


(C) 

Figure  4.9:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of  Plastic 
Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Plastic  Strain  as  Removal 

Criterion  (a/W  =  0.5) 
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step34  inci,  l=2.23e+001 
Plaslic  Stiain  Equivalent 


>  1.50e  +  00 

<  1.50e  +  00 

<  1.38e  +  00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-01 

<  1.15e-01 

<  O.OBe  +  OB 


Max  =  1.84e^BB 
Min  =  B.BBe  +  BB 


(a) 


(b) 


step34  incl,  l=2.23e+001 
Plastic  Strain  Equivalent 


>  1 .5Be  +  BB 

<  1.5Be+BB 

<  1.38e+BB 

<  1.27e+BB 

<  1.15e  +  BB 

<  1  .B4e  +  BB 

<  9.23e-Bl 

<  8.B8e-B1 

<  6.92e-B1 

<  5.77e-Bl 

<  4.62e-B1 

<  3.4Be-B1 

<  2.31  e-BI 

<  1.15e-B1 

<  B.BBe  +  BB 


Max  =  1.84e  +  BB 
Min  =  B.BBe  +  BB 


A 


z 


(c) 


Figure  4.10:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  3  Material  Properties  Using  Plastic  Strain  as 

Removal  Criterion  (a/W  =  0.4) 
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step50  inc1,  t=3.23e+001 
Plastic  Strain  Equivalent 


>  1 .50e  +  00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-01 

<  1.15e-01 

<  0.06e  +  08 


Max  =  1.92e  +  B0 
Min  =  0.08e  +  00 


A 


O- 

z 


(a) 


(b) 


stepSO  inci,  l=3.23e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-B1 

<  5.77e-81 

<  4.62e-01 

<  3.46e-B1 

<  2.31e-Bl 

<  1.15e-B1 

<  O.OOe  +  00 

Max  =  1.92e  +  0B 
Min  =  8.88e  +  B8 


A 


<3- 

z 


(c) 


Figure  4.1 1:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  3  Material  Properties  Using  Plastic  Strain  as 

Removal  Criterion  (a/W  =  0.45) 
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step69  inci,  l=4.35e+001 
Plastic  Strain  Equivalent 


>  1 .50e  +  00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-B1 

<  5.77e-B1 

<  4.62e-B1 

<  3.46e-B1 

<  2.31e-B1 

<  1.15e-B1 

<  8.08e  +  88 


Max  =  1.92e  +  BB 
Min  =  8.8Be  +  B8 


A 


(a) 


(b) 


step69  incl.  t=4.35e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-B1 

<  4.62e-01 

<  3.46e-81 

<  2.31e-01 

<  1.15e-B1 

<  O.OOe  +  00 

Max  =  1.92e  +  00 
Min  =  0.88e  +  08 


A 


<3- 

z 


(C) 


Figure  4.12:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  3  Material  Properties  Using  Plastic  Strain  as 

Removal  Criterion  (aAV  =  0.5) 
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step55  inci,  l=3.24e+001 
Plastic  Stiain  Equivalent 


>  1 .50e  +  00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-B1 

<  5.77e-B1 

<  4.62e-B1 

<  3.46e-B1 

<  2.31e-B1 

<  1.15e-B1 

<  0.08e  +  00 


Max  =  2.25e  +  B0 
Min  =  O.OOe  +  00 


A 


(a) 


(b) 


step55  inci,  t=3.24e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-01 

<  1.15e-01 

<  O.OOe  +  00 


Max  =  2.25e  +  00 
Min  =  O.OOe  +  00 


A 


C- 

z 


(C) 


Figure  4.13:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  1  Material  Properties  Using  Strain  Energy 

Density  as  Removal  Criterion  (a/W  =  0.4) 
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step95  inci,  t=5.15e+001 
Plastic  Strain  Equivalent 


>  1 .50e  +  00 

<  1.50e+00 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-Bl 

<  5.77e-81 

<  4.62e-81 

<  3.46e-01 

<  2.31e-01 

<  1.15e-81 

<  8.88e  +  88 


Max  =  2.25e  +  B0 
Min  =  O.OOe  +  00 


A 


O- 

z 


(a) 


(b) 


step95  inc1.  t=5.15e+001 
Plastic  Strain  Equivalent 


>  1.50e+00 

<  1.50e+80 

<  1.38e+00 

<  1.27e+00 

<  1.15e  +  00 

<  1 .04e  +  00 

<  9.23e-01 

<  8.08e-01 

<  6.92e-01 

<  5.77e-01 

<  4.62e-01 

<  3.46e-01 

<  2.31e-01 

<  1.15e-01 

<  O.OOe  +  00 

Max  =  2.25e  +  00 
Min  =  O.OOe  +  00 


(c) 


Figure  4.14:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  1  Material  Properties  Using  Strain  Energy 

Density  as  Removal  Criterion  (a/W  =  0.43) 
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step34  inci,  l=3.03e+001 
Plastic  Strain  Equivalent  (N) 


>  5.00e-01 

<  B.OOe-OI 

<  4.62e-01 

<  4.23e-01 

<  3.85e-01 

<  3.46e-01 

<  3.08e-01 

<  2.7Be-01 

<  2.31e-01 

<  1.93e-01 

<  1.55e-01 

<  1.16e-01 

<  7.79e-02 

<  3.95e-02 

<  1.11e-03 


Max  =  1.65e  +  00 
Min  =  1.11e-B3 


A 


C3- 

z 


(a) 


(b) 


slep34  inci,  l=3.03e+001 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  5.00e-01 

<  4.62e-01 

<  4.23e-01 

<  3.85e-01 

<  3.46e-01 

<  3.08e-01 

<  2.70e-01 

<  2.31e-01 

<  1.93e-01 

<  1.55e-01 

<  1.16e-01 

<  7.79e-02 

<  3.95e-02 

<  1.11e-03 


Max  =  1 .65e  +  00 
Min  =  1.11e-03 


A 


(c) 


Figure  4.15:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Strain  Energy 

Density  as  Removal  Criterion  (a/W  =  0.4) 
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step56  inci,  l=5.13e+001 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  B.OOe-OI 

<  4.62e-01 

<  4.24e-01 

<  3.B6e-01 

<  3.48e-01 

<  3.10e-01 

<  2.72e-01 

<  2.34e-01 

<  1.96e-01 

<  1.57e-01 

<  1.19e-01 

<  8.14e-02 

<  4.33e-02 

<  5.25e-03 


Max  =  1.39e  +  00 
Min  =  5.25e-03 


A 


(a) 


(b) 


step56  inci,  t=5.13e+001 
Plastic  Strain  Equivalent  (N) 


>  B.OOe-OI 

<  B.OOe-01 

<  4.62e-01 

<  A.2Ae-0^ 

<  3.86e-01 

<  3.48e-01 

<  3.10e-01 

<  2.72e-01 

<  2.34e-01 

<  1.96e-01 

<  1.B7e-01 

<  1.19e-01 

<  8.14e-02 

<  4.33e-02 

<  B.2Be-03 


Max  =  1.39e  +  00 
Min  =  B.2Be-03 


A 


O- 

z 


(C) 


Figure  4.16:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Strain  Energy 

Density  as  Removal  Criterion  (a/W  =  0.45) 
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Figure  4.17:  (a)  Plastic  Strain  Contours,  (b)  Crack  Propagation  and  (c)  Close-up  View  of 
Plastic  Strain  Contours  for  Model  With  Curve  2  Material  Properties  Using  Strain  Energy 

Density  as  Removal  Criterion  (a/W  =  0.5) 
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5. 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Summary  and  Conclusions 

This  study  focused  on  the  development  of  shear-lip  in  DT  specimens  under  quasi-static 
loading  rate.  The  objective  of  the  study  was  to  demonstrate  the  feasibility  of  a  numerical 
simulation  of  shear-lips  for  a  25mm  base  metal  DT  specimen  in  three-point  bend  test.  Three  real 
stress-real  strain  curves  representing  three  different  materials  with  varying  degrees  of  plastic 
behaviour  were  considered.  These  material  curves  are  designated  as  Curves  1,  2  and  3.  Curve  1 
represents  the  actual  real-stress-  real-strain  curve  of  the  350WT  steel  material.  On  the  other  hand, 
Curve  2  represented  the  same  curve,  but  with  the  fracture  strain  reduced  to  have  that  of  Curve  1 ; 
while  Curve  3  represents  a  similar  curve  as  Curve  1,  but  with  the  yield  stress  reduced  to  half  the 
Curve  1  value.  Crack  advance  and  element  or  node  release  was  based  on  two  possible  criteria, 
namely  (a)  when  the  average  strain  in  the  crack  tip  elements  that  occupy  a  physical  surface  of  1 .25 
mm  (0.050  inches)  exceeds  the  strain  at  fracture  on  the  supplied  data;  and  (b)  when  the  strain 
energy  density  in  elements  occupying  a  volume  defined  by  1 .25  mm  (0.050  inches)  of  free  surface 
at  the  crack  tip  and  1 .25  mm  (0.050  inches)  of  material  depth  exceeds  the  strain  energy  density 
calculated  from  the  supplied  real  stress  real  strain  curves. 

It  was  shown  that  the  plastic  strain  contour  patterns  obtained  from  the  finite  element 
analyses  could  be  used  to  estimate  the  shear-lip  sizes.  At  any  given  crack  size  (a/W  level),  the 
shear- lip  size  is  smallest  for  Curve  2  material  and  largest  for  Curve  3.  This  trend  is  reasonable  as 
the  Curve  2  material  is  the  least  plastic  (100%  strain  to  failure  as  opposed  to  200%  strain  for 
Curves  1  and  3);  and  the  Curve  3  material  is  the  most  plastic  (has  lower  yield  stress  of  197MPa 
compared  to  395MPa  for  Curves  1  and  2,  and  high  failure  strain  of 200%).  For  a  given  material, 
the  shear-lip  size  tends  to  reduce  with  increasing  crack  size.  For  instance,  for  the  Curve  1 
material,  the  shear-lip  sizes  for  a/W  =  0.4,  0.45  and  0.5  are  10.95mm,  10.30mm  and  9.70mm, 
respectively.  The  same  trends  of  shear-lip  sizes  were  observed,  regardless  of  whether  the  plastic 
strain  criterion  or  strain  energy  density  criterion  was  used  for  crack  advance.  However,  for  any 
given  material  and  crack  size,  the  shear-lip  size  predicted  based  on  the  strain  energy  density 
criterion  was  greater  than  that  predicted  based  on  the  plastic  strain  criterion. 
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5.2 


Recommendations 


It  is  recommended  that  comparisons  be  made  between  the  results  obtained  in  this  study 
and  experimental  studies  being  conducted  at  DRDC.  Based  on  these  comparisons,  it  might  be 
necessary  to  refine  the  methodology  presented.  In  this  regard,  issues  that  need  consideration 
include  investigation  of  the  influence  of  mesh  refinement  and  the  use  of  strain  energy  density 
criterion  that  is  obtained  from  tests  in  a  triaxial  state  of  stress.  It  is  also  recommended  that  the 
analysis  be  extended  to  specimens  of  different  thickness  (say  16mm  and  12.5mm).  Such  studies 
could  be  used  to  correlate  results  from  earlier  DRDC  studies  [4]  and  hence  lead  to  the  use  of  the 
shear-lip  size  information  as  a  material  property. 
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